Application of Boyd’s periodization and relaxation method in a spectral atmospheric limited-area model, part II : accuracy analysis and detailed study of the operational impact

Spectral limited-area models face a particular challenge at their lateral boundaries: the fields need to be made periodic. Boyd proposed a windowing-based method to improve the periodization and relaxation. In a companion paper, the implementation of this windowing method in the operational semi-implicit semi-Lagrangian spectral HARMONIE system was described and some first reproducibility tests, comparing this method to the old existing one, were presented. The present paper provides an in-depth study of the impact of this method for different configurations of the implementation. This is carried out in three steps in well-controlled experimental setups of increasing complexity. First, different aspects of Boyd’s method are analyzed in an idealized perfect-model test using a representative 1D shallow-water model. Second, the implementation is tested in an adiabatic 3D numerical weather prediction (NWP) model with perfect-model experiments. Finally, the impact of using Boyd’s method in a more operational-like NWP context is investigated as well. The presented tests show that, while the implementation of Boyd’s method is neutral in terms of scores, it is superior to the existing spline method in the case of strong dynamical forcings at the lateral boundaries

Simulating model uncertainty of subgrid-scale processes by sampling model errors at convective scales

Ideally, perturbation schemes in ensemble forecasts should be based on the statistical properties of the model errors. Often, however, the statistical properties of these model errors are unknown. In practice, the perturbations are pragmatically modelled and tuned to maximize the skill of the ensemble forecast. In this paper a general methodology is developed to diagnose the model error, linked to a specific physical process, based on a comparison between a target and a reference model. Here, the reference model is a configuration of the ALADIN (Aire Limitée Adaptation Dynamique Développement International) model with a parameterization of deep convection. This configuration is also run with the deep-convection parameterization scheme switched off, degrading the forecast skill. The model error is then defined as the difference of the energy and mass fluxes between the reference model with scale-aware deep-convection parameterization and the target model without deep-convection parameterization. In the second part of the paper, the diagnosed model-error characteristics are used to stochastically perturb the fluxes of the target model by sampling the model errors from a training period in such a way that the distribution and the vertical and multivariate correlation within a grid column are preserved. By perturbing the fluxes it is guaranteed that the total mass, heat and momentum are conserved. The tests, performed over the period 11–20 April 2009, show that the ensemble system with the stochastic flux perturbations combined with the initial condition perturbations not only outperforms the target ensemble, where deep convection is not parameterized, but for many variables it even performs better than the reference ensemble (with scale-aware deep-convection scheme). The introduction of the stochastic flux perturbations reduces the small-scale erroneous spread while increasing the overall spread, leading to a more skillful ensemble. The impact is largest in the upper troposphere with substantial improvements compared to other state-of-the-art stochastic perturbation schemes. At lower levels the improvements are smaller or neutral, except for temperature where the forecast skill is degraded

Jean-François Geleyn, professeur à l'université de Gand

Jean-François Geleyn had a strong link with Belgium through the activities of the ALADIN consortium. When Ghent university started a postgraduate program on meteorology and numerical weather prediction in 2007 he became a guest professor and one of the main driving forces of the creation of this academic program. Later the curriculum of this program was extented to include climate modeling. Many experts at the Royal Meteorological Institute of Belgium followed his courses. Recently, the activities at Ghent university formed the basis for the creation of the Atmospheric Physics Unit within the Department of Physics and Astronomy

Generalization and application of the flux-conservative thermodynamic equations in the AROME model of the ALADIN system

General yet compact equations are presented to express the thermodynamic impact of physical parameterizations in a NWP or climate model. By expressing the equations in a flux-conservative formulation, the conservation of mass and energy by the physics parameterizations is a built-in feature of the system. Moreover, the centralization of all thermodynamic calculations guarantees a consistent thermodynamical treatment of the different processes. The generality of this physics-dynamics interface is illustrated by applying it in the AROME NWP model. The physics-dynamics interface of this model currently makes some approximations, which typically consist of neglecting some terms in the total energy budget, such as the transport of heat by falling precipitation, or the effect of diffusive moisture transport. Although these terms are usually quite small, omitting them from the energy budget breaks the constraint of energy conservation. The presented set of equations provides the opportunity to get rid of these approximations, in order to arrive at a consistent and energy-conservative model. A verification in an operational setting shows that the impact on monthly-averaged, domain-wide meteorological scores is quite neutral. However, under specific circumstances, the supposedly small terms may turn out not to be entirely negligible. A detailed study of a case with heavy precipitation shows that the heat transport by precipitation contributes to the formation of a region of relatively cold air near the surface, the so-called cold pool. Given the importance of this cold pool mechanism in the life cycle of convective events, it is advisable not to neglect phenomena that may enhance it

Predicting small-scale, short-lived downbursts : case study with the NWP limited-area ALARO model for the Pukkelpop thunderstorm

The authors consider a thunderstorm event in 2011 during a music festival in Belgium that produced a short-lived downburst of a diameter of less than 100 m. This is far too small to be resolved by the kilometric resolutions of today's operational numerical weather prediction models. Operational forecast models will not run at hectometric resolutions in the foreseeable future. The storm caused five casualties and raised strong societal questions regarding the predictability of such a traumatic weather event.In this paper it is investigated whether the downdrafts of a parameterization scheme of deep convection can be used as proxies for the unresolved downbursts. To this end the operational model ALARO [a version of the Action de Recherche Petite Echelle Grande Echelle-Aire Limitee Adaptation Dynamique Developpement International (ARPEGE-ALADIN) operational limited area model with a revised and modular structure of the physical parameterizations] of the Royal Meteorological Institute of Belgium is used. While the model in its operational configuration at the time of the event did not give a clear hint of a downburst event, it has been found that (i) the use of unsaturated downdrafts and (ii) some adaptations of the features of this downdraft parameterization scheme, specifically the sensitivity to the entrainment and friction, can make the downdrafts sensitive enough to the surrounding resolved-scale conditions to make them useful as indicators of the possibility of such downbursts

Application of Boyd’s periodization and relaxation method in a spectral atmospheric limited-area model, part I: implementation and reproducibility tests

This paper describes the implementation of a proposal of Boyd for the periodization and relaxation of the fields in a full three-dimensional spectral semi-implicit semi-Lagrangian limited-area model structure of an atmospheric modeling system called HARMONIE that is used for numerical weather prediction and regional climate studies. Some first feasibility tests in an operational numerical weather prediction context are presented. They show that, in terms of standard operational forecast scores, Boyd’s windowing-based method provides comparable performance as the old existing spline-based periodization procedure. However, the real improvements of this method should be expected in specific cases of strong dynamical forcings at the lateral boundaries. An extensive demonstration of the superiority of this windowing-based method is provided in an accompanying paper

Evaluating the performance of SURFEXv5 as a new land surface scheme for the ALADINcy36 and ALARO-0 models

The newly developed land surface scheme SURFEX (SURFace EXternalisee) is implemented into a limited-area numerical weather prediction model running operationally in a number of countries of the ALADIN and HIRLAM consortia. The primary question addressed is the ability of SURFEX to be used as a new land surface scheme and thus assessing its potential use in an operational configuration instead of the original ISBA (Interactions between Soil, Biosphere, and Atmosphere) scheme. The results show that the introduction of SURFEX either shows improvement for or has a neutral impact on the 2m temperature, 2m relative humidity and 10m wind. However, it seems that SURFEX has a tendency to produce higher maximum temperatures at high-elevation stations during winter daytime, which degrades the 2m temperature scores. In addition, surface radiative and energy fluxes improve compared to observations from the Cabauw tower. The results also show that promising improvements with a demonstrated positive impact on the forecast performance are achieved by introducing the town energy balance (TEB) scheme. It was found that the use of SURFEX has a neutral impact on the precipitation scores. However, the implementation of TEB within SURFEX for a high-resolution run tends to cause rainfall to be locally concentrated, and the total accumulated precipitation obviously decreases during the summer. One of the novel features developed in SURFEX is the availability of a more advanced surface data assimilation using the extended Kalman filter. The results over Belgium show that the forecast scores are similar between the extended Kalman filter and the classical optimal interpolation scheme. Finally, concerning the vertical scores, the introduction of SURFEX either shows improvement for or has a neutral impact in the free atmosphere

The urban climate of Ghent, Belgium : a case study combining a high-accuracy monitoring network with numerical simulations

As urban environments have a specific climate that poses extra challenges (e.g. increased heat stress during heat waves), gaining detailed insight into the urban climate is important. This paper presents the high-accuracy MOCCA (MOnitoring the City's Climate and Atmosphere) network, which is monitoring the urban climate of the city of Ghent since July 2016. The study illustrates the complementarity between modelling and observing the urban climate. Two different modelling approaches are used: 1 km resolution runs of the SURFEX land surface model and 100 m resolution runs of the computationally cheaper UrbClim boundary layer model. On the one hand, urban models are able to simulate the spatial variability of the urban climate. As such, these models serve as a tool to help deciding on the locations of the measurement stations. On the other hand, the MOCCA observations are used to validate the high-resolution urban model experiments for the summer (July-August-September) of 2016. Our results demonstrate that the models capture the nighttime intra-urban temperature differences, but they are not able to reproduce the observed daytime temperature differences which are determined by the micro-scale environment